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Goal 

• Increase flexibility of T-MATS 

• Cantera increases flexibility of 
thermodynamics 

• Can model any flow 

• M-file elements allow users to prototype 
engineering elements 

• Slower than standard T-MATS 
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T-MATS 



• Simulinkcode 

• Library of thermodynamic elements 

- Standard library includes elements typical of aeropropulsion 

• Newton Raphson solver 

• Default thermodynamic table is air, water, and a 
hydrocarbon fuel 

• Systems can be modeled outside the standard 
elements/thermo 

- Create new thermo tables 

- Create elements 
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Cantera 



• Object-oriented software tools for problems involving 
chemical kinetics, thermodynamics, and/or transport 
properties 

• C++ based code with interfaces for python, matlab, 
C, and fortran 90 

• https://code.google.eom/p/cantera/ 
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Integration of T-MATS with Cantera 



• Allows any fluid combination to be modeled 

• Specify the thermodynamics of the possible products 

- Similarto CEA thermo. inp file 

• Requires specification of all “reactants” for the 
simulation 

- Similarto CEA reactant cards 

- Specify the different possible starting flows by composition 
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Species = { .7547 .232.0128 0 0 0; 


1 0000 0 ; 

.922189 .077811 0 0 0 0 0; 

00000 0 ; 00000 0 ; 00000 0 }; 
Name = { 'N2' '02' 'AR ; 


H20' 


ii ii ii ii ii 


l CH2"CH' 


ii ii ii ii 


ii ii ii ii ii a ii ii ii ii ii a ii ii ii ii ii ii 
3 3 

• Species and Name arrays need to be defined 


• A model with this definition can run with mixtures of Air, 
Water, and JP-7 


• Allows for models of aircraft engines with humidity 
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T-MATS Cantera Fluid Arrays 



Information 

Index 

Description 

W 

1 

Weight of the flow 

Tt 

2 

Total temperature 

Pt 

3 

Total pres sure 

ht 

4 

Total enthalpy 

comp 1 (to) comp 10 

5-14 

Percentageofflowcompositionforreactants lto 10 

s 

15 

Entropy 

rhot 

16 

Total density 

Ts 

17 

Static temperature 

Ps 

18 

Static pressure 

hs 

19 

Static enthalpy 

rhos 

20 

Static density 

Vflow 

21 

Flow velocity 

MN 

22 

How Mach number 

A 

23 

How area 

gamt 

24 

Totalgamma 

gams 

25 

Static gamma 


Each fluid location in a thermodynamic model is 
represented by an array that contains all the fluid 
properties at a given location 


www.nasa.gov 


National Aeronautics and Space Administration 


f-MATS Cantera Fluid Functions 



Function 

Description 

add ( flowl , flow2) 

Add add flowl and flow2 together, conserving enthalpy and 
mass 

copyFlow ( flow) 

Copy the information fromj How to anotherflow 

getMass Fraction ( flow, c ) 

Return the mass fraction of compound c intheobject^Zow 

set hP ( flow, ht , Pt ) 

Set the total conditions based on flow , total enthalpy and total 
pressure 

set MN1 (flow) 

Set the static conditions to sonic based on flow conditions 

set MNP s (flow, Ps) 

Set the static conditions based on flow and input static 
pressure 

set SP ( flow, S , Pt ) 

Set the total conditions based on flow and input entropy and 
total pres sure 

set TP ( flow, Tt , Pt ) 

Set the total conditions based on flow, total temperature and 
total pres sure 

set Ts P sMN ( fl o w, T s , P s , MN ) 

Set the conditions based on flow, static temperature, static 
pressure andMach 


• All communication between Cantera and T-MATS is 
handled by these functions 

• Functions return a new Cantera Fluid Array based on 
inputs (see previous slide) 
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f-MATS Element Files 



• Library of standard elements released in Simulink m- 
file format 

• Allows for development and prototyping 

• Elements are interpreted 

- No need to compile 

• Engineers can quickly create new elements 

• Block sets are released with T-MATS Cantera 
package 
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Instance Information 

• Needed a way to store instance information from one 
pass to another 

• Created to functions to store and retrieve information 
from one pass to another 

- Variables are stored in the MATLAB workspace with the 
object instance name attached to the variable instance name 

• setV sets the value or a variable in the workspace 

• getV gets the value of a variable from the workspace 

path = stripchar ( gcb ( ) ); 

setV ( 's_C_Nc', path, s_C_Nc ); 

s C Nc = getV ( ' s C Nc', path ); 
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Some Examples from Compressor Element 

• Setting the exit conditions 

FOideal = set_SP ( FI, FI(s),PtOut ); 
htOut = FI (ht) + ( FOideal (ht) - FI (ht) )/eff 

% set the exit conditions to known enthalpy and 
%pressure 

FO = set hP ( FI, htOut, PtOut ); 
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Some Examples from Compressor Element 



• Design Point Scaling 


if IDes < .5 

s_eff = effDes / effMap; 
s_PR = ( PRdes - 1 )/( PRmap - 1 ); 
s_Wc = Wcln/ WcMap; 
setV( 's_eff, path, s_eff ); 
setV( 's_Wc', path, s_Wc ); 
setV( 's_PR', path, s_PR ); 
elseif IDes < 1 .5 

% get the maps scalars from the workspace 
s_eff= getV( 's_eff, path ); 
s_Wc= getV( 's_Wc', path ); 
s_PR= getV( 's_PR', path ); 
else 

% use the input values 
s_eff = s_eff Jn; 
s_Wc = sWcJn; 
s_PR = s_PR_in; 
end 
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Simulink Objects 




Mask appearance 

Describes port labels and colors 

Label colors are standard based on T-MATS style 
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• Parameter list 

• Lists the variables that can be input by the 
user to the dialog box 
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• S-function block parameters 

• Utilizes m-file to create S-function 

• Maps parameter dialog box to m-file 
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furbofan Model -JT9D 
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furbofan Model -JT9D 




NPSS with JANAF 
Output 

TMATS Cantera 
Output 

Altitude 

34000 ft 

34000 ft 

Mach 

number 

.8 

.8 

Weight flow 

674 lbm/sec 

674 lbm/sec 

Thrust 

11194 lbf 

11182 lbf 

SFC 

.6113 

.6116 
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Fuel Cell Model 



I tSOFC 




File Edit 

View Display Diagram Simulation Analysis Code Tools Help 


^ & 



▼ 

K ▼ 10.0 [Normal 

- 


Model Browser 



tSOFC 


© tSOFC 


Ospte,3 



100% 


ode45 


• Reactants are Air, H2, 02, and H20 
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Fuel Cell Model 



Specifying the reactants: 


Species ={ .7547 .232.0128 0 0 0; 
1 0 0000 : 


100000 

100000 


00000 0;00000 0 }; 


Name = { 'N2' '02' 'AR ; 

!|_I2 | n m n n m > j Q2' •• 

M n 11 11 n ii ii 


ii ii ii ii. 

5 

II II K, II II II II II II 

5 


Getting the mass fractions of an element: 

xN2_cOut= getMassFraction( FI 02, 'N2' ); 
x02_c0ut = getMassFraction( FI 02, '02' ); 
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Fuel Cell Model 

Removing oxygen from the flow: 

x02_Cathode1 = getMassFraction( FI Cathodel , '02') 

%Composition as mass flow (g/sec) 
w02_Cathode 1 = x02_Cathode 1 * w_Cathode 1 

%Composition as molar flow rate (mol/sec) 
M_02_Cathode1 = w02_Cathode1 / 32. 



%Calculates composition after electrochemistry... 

M_02_Cathode2 = M_02_Cathode1 - ((M_H2_Anode1 / 2.0) * pctH2util); 

M_02_Cathode2 = M_02_Cathode1 - ((M_H2_Anode1 / 2.0) * pctH2util); 
w02_lost = (M_02_Cathode1 - M_02_Cathode2)*32. * 0.002205 % Ib/sec 

Fl_temp02(8) = 1 ; 

FI_temp02(W) = -w02_lost; 

Fl_temp02= set_TP( Fl_temp02, FI_Cathode2(Tt), FI_Cathode2(Pt) ); 
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Conclusion 

• Cantera has been integrated with T-MATS 

• Capable of modeling any thermodynamic flow 

• Simulink block sets and MATLAB m-files 

• Allows for prototyping 

• Greatly increases the flexibility of T-MATS 

• Slower than standard T-MATS 
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• Download information may be found at: 
https://github.com/nasa/T-MATS/releases/ 
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